Magnetic ordering in the static intermediate-valent cerium compound Ce 2 RuZn 4 
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The low-temperature behavior of Ce2RuZn4 has been investigated. Specific heat and magnetic 
susceptibility data reveal an antiferromagnetic transition at a Neel temperature of 2 K. Ce2RuZn4 
is a static intermediate-valent compound with two crystallographically independent cerium atoms. 
The magnetic data clearly show that only one cerium site is magnetic (Ce 3+ ), while the second one 
carries no magnetic moment. The experimental data are interpreted with the help of first principles 
electronic structure calculations using density functional theory and the augmented spherical wave 
method. The calculations reveal the occurrence of two different cerium sites, which are characterized 
by strongly localized magnetic moments and strong Ce-Ru bonding, respectively. 

PACS numbers: 61.10.Nz, 75.50.Cc, 75.30.Sg 
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I. INTRODUCTION 

Ternary intcrmetallic compounds in the system Ce-Ru- 
X, where X is an element of the 3rd, 4th, or 5th main 
group, have been intensively studied in recent years. This 
is due to the fact that these materials exhibit peculiar 
structural and physical properties, where the transition 
metal ruthenium seems to play a key role. In compar- 
ison to related compounds with other transition met- 
als, ruthenium often shows extremely short Cc-Ru dis- 
tances, even much shorter than the sum of the cova- 
lent radii of 289 pmi Already the binary Laves phase 
CcRu2, 2 which is a superconductor below 6.2 is ex- 
traordinary In this compound, the cerium atoms are in 
an intermediate valent state and the a lattice parame- 
ter is much smaller than that of NdRti2 . When it comes 
to the ternary compounds, these effects are even more 
pronounced. The recently reported indium compounds 
CcRu .88lii2, Cci 6 Ru 8 In 3 7, Cc 3 Ru 2 In2, and Ce 3 Ru 2 In 3 
exhibit extremely short Ce-Ru distances^^ii which are 
directly associated with a strong tendency of some of the 
cerium sites towards tetravalency. 

The stannide CeRuSn also shows this peculiar 
behaviori^ The room-temperature modification adopts 
a superstructure of the CeCoAl type and the two crys- 
tallographically independent cerium sites are ordered on 
distinct positions. In particular, the intermediate valent 
Cel atoms have very short Cel-Ru distances of only 233- 
246 pm. Consequently, our electronic structure calcula- 
tions revealed strong Ce-Ru bonding. 

Another example for trivalcnt/intermcdiatc valent 
cerium ordering is the structure of Ce2RuZn4^ This 



structure type crystallizes in the space group PA/nmm 
with two formula units per cell and two twofold cerium 
sites with distinctly different coordination. Preliminary 
magnetic susceptibility investigations underlined the idea 
of trivalcnt/intermcdiatc valent cerium ordering. Ex- 
tending our systematic studies of such ruthenium based 
intermetallic compounds we here report on a compre- 
hensive investigation of the low-temperature properties 
of Ce2RuZii4. In particular, our work includes X-ray 
diffraction, specific heat and magnetic susceptibility mea- 
surements as well as electronic structure calculations. 



II. EXPERIMENTAL 

A. Synthesis 

Starting chemicals for the preparation of CC2RUZ114 
were a cerium ingot (Johnson Matthey, purity > 99.9 %), 
ruthenium powder (Degussa-Hiils, ca. 200 mesh, purity 
> 99.9 %), and zinc granules (Merck, purity > 99.9 %). 
Pieces of the cerium ingot were first arc-melted to a small 
button under argon The argon was purified before with 
molecular sieve, silica gel and titanium sponge (900 K). 
The elements were then mixed in the 2:1:4 atomic ratio 
(ruthenium powder was cold-pressed to a pellet of 6 mm 
diameter) and arc-welded in a tantalum ampoule under 
an argon pressure of 700mbar. The ampoule was sub- 
sequently placed in a water-cooled sample chamber of a 
high-frequency furnace (Huttinger Elektronik, Freiburg, 
Typ TIG 5/300) The heat treatment was similar to our 
previous experiments^ three times for lmin at 1360 K, 
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for lOmin at 1360 K, rapid cooling to 920 K, and a fi- 
nal annealing at that temperature for another 2 hours 
followed by quenching. 

B. Powder X-ray Diffraction 

The purity of the sample was checked through a 
Guinier powder pattern using Cu-K Q i radiation and a- 
quartz (a = 491.30 pm, c = 540.46 pm) as an internal 
standard. The tetragonal lattice parameters were ob- 
tained from a least-squares refinement of the Guinier 
powder data. The refined lattice parameters (a = 
719.6(1) pm, c = 520.0(1) pm) are in good agreement 
with the previously reported data (a = 719.6(1) pm, c = 
520.2(1) pm)^ No impurity phases have been detected 
on the level of X-ray powder diffraction. 

C. Specific Heat and Magnetic Susceptibility 

The temperature-dependent heat capacity was inves- 
tigated by means of a quasi-adiabatic step heating tech- 
nique employing a Quantum Design PPMS system. The 
measurements were performed on a polycrystalline sam- 
ple thermally connected by grease (Apiezon-N) at tem- 
peratures ranging from 1.8 up to 300 K. For temperatures 
below 2K, down to about 80mK, specific heat data were 
also collected in a 3 He/ 4 He-dilution-cryostat using a re- 
laxation method The uncertainty of the measurements 
reported in this paper is estimated to be below 3%. The 
magnetic measurements in the temperature range 1.8 K 
< T < 300 K were done using a Quantum Design MPMS7 
superconducting quantum interference device (SQUID) 
magnetometer. 

III. THEORETICAL METHOD 

The calculations are based on density-functional the- 
ory (DFT) and the local density approximation (LDA). 
They were performed using the scalar-relativistic im- 
plementation of the aug mented spherical wave (ASW) 
method (see Refs. IT5lll6l and references therein). In the 
ASW method, the wave function is expanded in atom- 
centered augmented spherical waves, which are Hankel 
functions and numerical solutions of Schrodinger's equa- 
tion, respectively, outside and inside the so-called aug- 
mentation spheres. In order to optimize the basis set, 
additional augmented spherical waves were placed at 
carefully selected interstitial sites. The choice of these 
sites as well as the augmentation radii were automati- 
cally determined using the sphere-geometry optimization 
algorithm^ Self consistency was achieved by a highly ef- 
ficient algorithm for convergence acceleration.— The Bril- 
louin zone integrations were performed using the linear 
tetrahedron method with up to 105 k-points within the 
irreducible wedge i 16 i 19 



In the present work, we used a new full potential ver- 
sion of the ASW method, which was implemented only 
very recently^ In this version, the electron density and 
related quantities are given by a spherical harmonics ex- 
pansion inside the muffin-tin spheres. In the remaining 
interstitial region, a representation in terms of atom- 
centered Hankel functions is uscd*2i However, in contrast 
to previous related implementations, we here get away 
without needing a so-called multiple-K basis set, which 
fact allows for a very high computational speed of the 
resulting scheme. 

IV. RESULTS AND DISCUSSION 

A. Crystal Chemistry 

Cc2RuZn4 crystallizes with a peculiar structure type 
with space group PA/nmm^ Since the crystal chem- 
istry of this zinc-rich intermetallic compound has been 
described in detail in our previous work^ here we fo- 
cus only on the cerium coordination, which is responsible 
for the distinctly different magnetic behaviour of the two 
cerium sites. A view of the Ce2RuZn4 structure approx- 
imately along the c axis is presented in Fig. [T] The zinc 
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FIG. 1: The crystal structure of Ce2RuZri4: (left-hand side) 
view of the structure approximately along the c axis and 
(right-hand side) coordination polyhedra and site symmetries 
of the two crystallographically independent cerium sites. The 
cerium, ruthenium, and zinc atoms are drawn as medium grey, 
filled, and open circles, respectively. 

atoms build up a complex three-dimensional network, 
which leaves two different kinds of channels parallel to 
the c axis. They are filled by the cerium and ruthenium 
atoms. As emphasized in the right-hand part of Fig. [TJ 
the trivalent Cel atoms have 12 nearest neighbours of 
only the zinc type. In contrast, the intcrmediate-valent 
(nearly tetravalcnt) Ce2 atoms have a smaller coordina- 
tion number of 10 atoms with two ruthenium and eight 
zinc sites. However, the most striking structural feature 
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are the short Ce2-Ru distances of 260 pm, which are even 
shorter than the sum of the covalent radii of 289pniji 
This seems to be a characteristic of ruthenium and leads 
to a destabilization of the cerium valence. The peculiar 
electronic properties, which arise from this unusual crys- 
tal chemistry, are addressed below. 



B. Specific Heat 

In Fig. [21 the molar specific heat capacity divided 



Subtracting the calculated electronic, lattice and nu- 
clear parts from the total specific heat data results in the 
magnetic contribution of specific heat C mag . Subsequent 
integration of C mag /T leads to the magnetic entropy ac- 
cording to: 

Smag = / ^p<dT = R\n(2S + 1) . 

Here, S mag and S refer to the magnetic entropy and the 
spin, respectively. The corresponding graph is shown in 
the inset of Fig. [2] The resulting magnetic entropy of 
i?ln2 is due to a two level 5 = 1/2 state. 
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FIG. 2: (Color online) Specific heat divided by tempera- 
ture vs. logT of Ce2RuZn4 in the temperature range between 
50 mK and 300 K. The slight increase of C/T below IK is 
mainly due to the nuclear quadrupole moments displayed by 
the 101 Ru and 67 Zn isotopes; the solid red line represents the 
lattice and electronic contributions. The inset shows the mag- 
netic entropy. The red line displays Smag — -Rln2, which is 
the limit for a two level system. 



C. Magnetic Susceptibility 

The magnetic configuration of the Ce ion can be 
proven by susceptibility measurements O ur sus- 

ccptibility data collected at a field of 0.5 T are shown 
in Fig. [3] Here the molar core diamagnctism of Xdia = 
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by temperature C/T of Ce2RuZn4 is displayed in a 
semi-logarithmic plot for temperatures between 0.08 and 
300 K. A clear peak at 2K indicates a pronounced mag- 
netic contribution to the heat capacity. In order to ex- 
tract the magnetic contribution of the specific heat we 
determined the lattice contributions theoretically using 
a single Debye term and two Einstein modes. In this cal- 
culation, we fixed the number of internal degrees of free- 
dom to 21, according to the seven atoms in the unit cell. 
Using a weight distribution of 1:2:4, we obtained a De- 
bye temperature of Qd = 96 K, and two Einstein modes 
of 6 E i = 129 K and 9 E 2 = 238 K. The solid red line in 
Fig- El shows the corresponding result taking into account 
an electronic contribution of 7 = 30mJ/molK 2 . Below 
1 K, the specific heat of Ce2RuZn4 increases slightly giv- 
ing rise to a nuclear electronic Schottky contribution. 
We have calculated the hyperfinc contributions to the 
heat capacity originating in zero magnetic field mainly 
from the quadruple moments of 101 Ru and 67 Zn leading 
to an average local internal electrical field gradient of 
V zz = (5 ± 0.8) x 10 22 V/m 2 . 



FIG. 3: (Color online) The magnetic susceptibility x vs. logT 
of Ce2RuZn4. At low temperature the beginning of an antifer- 
romagnetic transition is observed. The inset exhibits the in- 
verse susceptibility 1/x vs. T. The solid line is a least squares 
fit according to the modified Curie- Weiss-law. 

—0.221 memu/mol is substracted from the data. Xdia is 
calculated using the expression Xdia = —0.79 x JT Zi x 
I0~ 6 (r /clb) 2 emu/mol, where JT Zi is the sum of all core 
electrons per formula unit, (r/as) is roughly estimated 
to be unity, and as is the Bohr radius. In the inset of Fig. 
[31 a linear temperature dependence of x _1 (T) is observed 
above 50 K. A fit according to the modified Curie Weiss- 
law 

C 



yields a temperature independent susceptibility xo = 
6.6 memu/mol, an effective moment /U e // = 2.57/i^ 
and the paramagnetic Curie- Weiss temperature Qcw = 
— 2.6K. The effective magnetic moment is close to the 
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theoretical value associated with the Ce 3+ state, i.e. 
2.54/is, but not to the calculated value of /z e // = 3.59/is 
as it is expected for the two Ce 3+ ions in the unit cell. 
These data clearly indicate that in Ce2RuZn4 the Ce ions 
occur in two different valence states Ce 3+ and Ce 4+ . A 
slight curvature of the magnetic susceptibility observed 
below 3 K is attributed to the beginning of an antiferro- 
magnctic transition with a Neel-temperature of Tn = 2 K 
as already derived from the specific heat data. 



D. Electronic Structure 

The investigation of the electronic structure was per- 
formed in two steps. While in a first set of calcula- 
tions, spin-degeneracy was enforced, we considered spin- 
polarization only at a later stage. This allowed to investi- 
gate the influence of the magnetic order on the electronic 
states and the total energy in more detail. 

The partial densities of states (DOS) resulting from 
the first step are displayed in Figs. 0] and According to 
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FIG. 4: (Color online) Partial densities of states of non- 
magnetic Ce2RuZn4. 

Fig. [31 the densities of states is dominated by the cerium 
4/ states. The different behaviour of the two cerium 
sites is clearly visible. While the 4/ states of the Cel 
atoms, which are at the centers of the Zn cages, form a 
very sharp peak right above the Fermi energy, the linearly 
coordinated Ce2 Af states give rise to rather broad bands 
of about 1 eV width due to their overlap with the Ru 5d 
states. 

The partial DOS due to the cerium and ruthenium 
d states are displayed in Fig. \5\ As expected, the Ru 
Ad states are located well below the Fermi energy. In 
contrast, the Ce 5d states lead to broad bands starting 
just slightly below Ep. Yet, hybridization of both types 
of bands leads to considerable contributions of the cerium 
and ruthenium states, respectively, below and above Ep. 
In particular, we observe a pronounced similarity of the 
Ru Ad partial DOS with that due to the Ce 5d states, 
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FIG. 5: (Color online) Partial densities of states of non- 
magnetic Ce2RuZn4. 



which is a signature of strong d-d overlap, which is even 
stronger for the Ce2 sites. In passing, we mention the 
Zn 3d states, which lead to a narrow band in the energy 
range from —8 to — 7eV and are not included in the plot. 

For simplicity reasons, subsequent spin-polarized cal- 
culations were performed for an assumed ferromagnetic 
order. These calculations resulted in a lowering of the to- 
tal energy by about 0.9 mRyd and a magnetic moment of 
about O.AfiB per formula unit, respectively Interestingly, 
the latter is carried by the Cel atoms alone, whereas the 
Ce2 states are not polarized. This becomes obvious from 
Fig. [51 which displays the partial densities of states of 
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FIG. 6: (Color online) Partial densities of states of assumed 
ferromagnetic Ce2RuZn4. 

assumed ferromagnetic Ce2RuZii4. Obviously, the finite 
spin-polarization at the Cel sites, which is mainly due 
to the 4/ states, results from the near- .Bp position of the 
latter. In contrast, the Ce2 4/ states are not polarized 
due to their much lower occupation. 

Motivated by the strongly localized nature of the Ce 4/ 
states, additional LDA+U calculations were performed 
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with U = 6.8 eV and J = 0.68 eV applied at both cerium 
sites. These calculations resulted in a magnetic moment 
of I.O^b at each Cel site as is expected for Ce 3+ . In 
contrast, the Ce2 sites remain unpolarized. In addi- 
tion, a complete restructuring of the electronic states 
as compared to the LDA calculations is observed in the 
partial densities of states shown in Fig. [71 While the 
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FIG. 7: (Color online) Partial densities of states of assumed 
ferromagnetic Ce2RuZn4. 

Ce2 4/ states remain centered at about 1 eV, the Cel 
4/ states experience a considerable splitting of about 
6eV due to the strong electronic correlations. While the 
spin-minority states are essentially empty, spin-majority 
states display a finite DOS at about — 2cV, which leads 
to the magnetic moments at this site. In contrast, the d 
states of all species remain essentially unpolarized. 

Finally, in order to investigate the differences in chem- 
ical bonding, we display in Figs. [5] and [5] the partial cova- 
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FIG. 8: (Color online) Partial covalence energies of assumed 
ferromagnetic Ce2RuZn4. 

lence energies E cov of assumed ferromagnetic Ce2RuZn 4 . 
In doing so, we concentrate on the bonding between the 



> 
c 
u 




-2 2 

(E - E F ) (eV) 



FIG. 9: (Color online) Partial covalence energies of assumed 
ferromagnetic Ce2RuZn4. 



Ru Ad states with the 5d and the 4/ states at the two 
different cerium sites. Note that negative and positive 
contributions to E cov point to bonding and antibonding 
states, respectively. The results give a very clear and 
impressive indication of the strong bonding between the 
Ru atoms and the Ce2 sites; the E cov curves are positive 
throughout up to Ep. In contrast, there is almost no 
bonding between the Ru atoms and the Cel sites, which 
are centered in the zinc cages. 



V. CONCLUSION 

The intcrmediate-valent compound Cc2RuZn4 was in- 
vestigated using specific heat and susceptibility measure- 
ments as well as first principles electronic structure cal- 
culations. This material is characterized by the presence 
of two different cerium sites with trivalent/intermediate 
valence ordering. According to the calculations, these dif- 
ferences go along with a striking variance in Ce-Ru bond- 
ing, which is strong at the intermcdiate-valent cerium 
sites and rather weak at the Cel positions at the cen- 
ter of the characteristic zinc cages. As a consequence, 
the LDA+U calculations reveal almost perfectly localized 
magnetic moments of 1.0/ib at the latter sites, whereas 
the Ru-bonded cerium displays no spin polarization. The 
theoretical findings arc in full agreement with the exper- 
imental data, which point to an antiferromagnetic tran- 
sition at a N eel-temperature of Tn = 2K. 
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